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BOX DIMENSION OF A HYPERBOLIC SADDLE LOOP 


MAJA RESMAN 


Abstract. We compute the box dimension of a spiral trajectory around a 
hyperbolic saddle loop, as the simplest example of a hyperbolic saddle poly cy¬ 
cle. In cases of weak foci and limit cycles, Zubrinic and Zupanovic show that 
the box dimension of a spiral trajectory is in a bijective correspondence with 
cyclicity of these sets. We show that, in saddle loop cases, the box dimension 
is related to the cyclicity, but the correspondence is not bijective. 

In addition, complex saddles are complexifications of weak foci points, as 
well as of hyperbolic saddles. Computing the box dimension around the sad¬ 
dle point of a hyperbolic saddle loop is hopefully a preliminary technique for 
computing the box dimension of leaves of a foliation around resonant complex 
saddles. 
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1. Introduction 

1.1. Motivation. Elementary limit periodic sets include weak foci, limit cycles 
and hyperbolic saddle polycycles. For a good overview of what is known about the 
cyclicities of these sets, see the book of Roussarie m- 

The overall use of fractal dimensions in dynamics is explained in short in e.g. |21j . 
One of the fractal dimensions is the box dimension. There is a relation between 
the box dimension and the cyclicity of limit periodic sets, noted by a group of 
authors in m , gsi, m • It was shown that the box dimension of either a spiral 
trajectory in the plane or of a discrete orbit of the Poincare map on the real line 
uniquely reveals the cyclicity of a limit cycle or a weak focus (in generic unfoldings). 
Box dimension of a trajectory is closely related to the asymptotic behavior of the 
Lebesgue measure of the £-neighborhood of the trajectory, as £ —> 0, see Section [T2| 
below for the precise definition. Only one (any) spiral trajectory is needed. 

In (2(7], box dimension of spiral trajectories around weak foci and limit cycles was 
computed. The box dimension was computed applying either the flow-box theorem 
PJ or the flow-sector theorem |20j , which state a locally parallel structure of the tra¬ 
jectory at a transversal to the set. Box dimension is computed as Cartesian product 
dimension, using the results about the box dimension of the one-dimensional orbit 
of the Poincare map given in [2]. The box dimension of a spiral trajectory is thus 
brought in a bijective correspondence with cyclicity in generic bifurcations. 

Our study here is concerned with relating geometric, fractal properties of one 
trajectory around the loop with the cyclicity of the loop. In this paper, we compute 
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the box dimension of a spiral trajectory around simplest hyperbolic saddle poly¬ 
cycles, hyperbolic loops, and compare it with cyclicity results. The descripton of 
bifurcations and the cyclicity of homoclinic sets, including hyperbolic saddle loops, 
is of interest in many applications, see the book of Palis, Takens [I2J. In JT5] or 
[7], the Poincare maps of hyperbolic saddle polycycles are described as having an 
expansion in a Chebyshev power-log scale. The cyclicity was related to the first 
non-zero coefficient in the expansion called the loop quantity , which is the analogon 
of Lyapunov coefficient for focus points. 

In cases of limit cycles and foci, the Poincare map has the same asymptotics on 
every transversal to the set. Here we show that the Poincare map at any transver¬ 
sal through the saddle vertex has a different asymptotics than Poincare maps at 
other transversals to the loop. Therefore, to compute box dimension of the spiral 
trajectory, we develop a new version of the flow-box theorem that identifies the 
trajectory around a saddle point with an appropriate set of hyperbolas. 

Furthermore, the Poincare map of limit cycles and foci is differentiable at zero. 
That is, their asymptotic expansion is a formal Taylor series. The Poincare maps 
of saddle polycycles expand in nondifferentiable real power-log Chebyshev scales, 
see IK], or Dl for Chebyshev scales. We show here that, unlike in differentiable 
cases from J2D], in nondifferentiable saddle-loop cases the box dimension is not 
sufficient for recognizing cyclicity. By its definition, the box dimension compares 
the area of the e-neighborhood to a power scale. It is implicitely shown in this paper 
that, in nondifferentiable cases, the areas of the e-neighborhoods of trajectories are 
not of power-type behavior. To get precise information, we need beforehand a 
finer logarithmic scale, depending on the bifurcation. The box dimension of a 
whole trajectory thus reflects the same deficiency as the box dimension of one- 
dimensional orbits of the Poincare map treated in HU: it cannot recognize between 
two neighboring cyclicities. 

We expect future applications of the result to be twofold. First, a similar tech¬ 
nique can be adapted to computing the box dimension of spiral trajectories around 
degenerate and nilpotent foci. As in the saddle loop case, the asymptotics of the 
Poincare map is dependent on the choice of the transversal. That is, there exist at 
least two distinct asymptotics. 

Secondly, planar fractal analysis obviously distinguishes between 
strong and weak foci, see m, m, but fails in distinguishing between hyperbolic 
saddle points. Namely, such points are not monodromic, trajectories are not re¬ 
curring and the Poincare map is not well-defined. Resonant complex saddles are 
complexifications of weak foci points, as well as of hyperbolic saddles. They are 
monodromic. Their monodromy maps are called holonomy maps , see e.g. m, 
[S]. The holonomy maps are complex parabolic germs defined on two-dimensional 
transversals. The box dimension of such germs was computed in [13] . By [16] , the 
structure is again locally parallel, except at the saddle vertex. Computing the box 
dimension of planar hyperbolic saddle loops is hopefully a preliminary technique 
for computing the box dimension of a leaf of a foliation around a complex saddle 
point. The box dimension is expected to reveal the formal type of the complex 
saddle. For our detailed conjecture on the box dimension of a complex resonant 
saddle and its relation to saddle formal invariants, see jUl Chapter 2]. 

Overview of the article. In Section [I] we state definitions and known results 
that we will use in the article. Section [2] is dedicated to the two main results, 


BOX DIMENSION, SADDLE LOOP 


3 


Theorem [l] and Theorem [2] Theorem [I] describes the asymptotic behaviour of the 
Poincare map on a transversal through the hyperbolic loop vertex. Theorem [2] 
gives the box dimension of a hyperbolic saddle loop. The results are proven in 
Section [3j In Section [0 we discuss two main applications. Subsection |d.l| explains 
the relation between the box dimension and cyclicity for a hyperbolic saddle loop. 
In Subsection |4.2[ we give a conjecture on the box dimension of a leaf of a foliation 
at a resonant complex saddle, related to the first formal invariant of the saddle. 


1.2. Definitions and notation. Let X be an analytic planar vector held with a 
monodromic hyperbolic saddle connection. That is, two eigenvalues of the linear 
part at the saddle point are of the opposite sign. After some rotation/translation 
changes of variables and possibly rescaling the time variable we can, without loss 
of generality, suppose that the loop is attracting, that the saddle lies at the origin 
and its separatrices correspond to the coordinate axes. The operations mentioned 
above do not change the phase portrait (only rotate and translate the loop), and 
therefore the box dimension of trajectories remains the same. Locally at the origin, 
the vector held is of the following form: 


x = x + P(x,y), 
y =-r-y + Q(x,y). 


Here, 0~ > 1 is the hyperbolicity ratio and P, Q are analytic functions of higher 
order than linear. 

If r £ Q?j_, the saddle is called resonant. We put r = p/q, where p, q £ N, 
( p , q) = 1. In the case r £ R+ \ Q, the saddle is nonresonant. 

Let T\ = {y = 1} and rj s {i = 1} represent the horizontal and the vertical 
transversal, parametrized so that the origin lies on the loop. Let P(s) denote the 
Poincare map on any transversal r not passing through the origin. The asymptotics 
(as s — y 0) of Poincare map is well-known and can be found in e.g. |151 Sections 
5.1.3, 5.2.20 


( 1 ) 


P(s) — s ~ s k or s k (— logs), k > 2, k £ N, if r = 1, 
P(s) ~ s r , if r > 1. 


The Poincare map s K > P{s) is a composition of the Dulac map s <—> D(s) around 
the saddle and the regular transition map around the loop s H > R(s). The leading 
term, if logarithmic s k (— logs), stems from the Dulac map, whereas the leading 
regular term s k stems from the regular loop transition map. 

The notion of codimension of the saddle loop is taken from [Hil Definition 27], 
as the number of the conditions imposed on the loop. In the case r = 1, the saddle 
loop is said to be of codimension 2k if s — P(s) ~ s fc , and of codimension 2k + 1 if 
s-P(s) ~ s fc+1 (— logs), for some k > 1. If r ^ 1, the saddle loop is of codimension 
1 . 


1 1n the case r < 1, the loop is repelling , so we rescale the time variable t o rl, and change 
the role of x and y. In this way, we get an attracting loop with r > 1 and the phase portrait 
remains the same. 

2 The notation ~ stands for: 

f (x} 

f ~ g, as x —»• 0, if lim —— = L. L / 0. 
x —.*-0 g(x) 
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In our computations in Section [2j we use the flow-box theorem from [2 p.75] or 

DO- 

The flow-box theorem. Let us consider a planar vector field of class C 1 . Assume 
that U C R 2 is a closed set the boundary of which is the union of two trajectories 
and two curves transversal to trajectories. If U is free of singularities and periodic 
orbits, then the vector field restricted to U is diffeomorphically conjugated to the 
field: 



on the unit square {(a;,?/)| 0 < x,y < 1}. That is, the flow on U can be represented 
as a parallel flow. 

We close the section with definitions of lower and upper box dimension and 
some properties used in the article. Let U C and its e-neighborhood, denoted 
U e , both be Lebesgue measurable. We denote the Lebesgue measure of the e- 
neighborhood by \U e \. The lower box dimension dim p (17) of U is defined as 

dim p (LO = infjs > 0 I liminf -0^- = 0}. 

Analogously, we define the upper box dimension of U, dims([/), with limsup in the 
above formula. If dim B (C/) = dim b{U), then 

chined) = dim P (U) = dims([/) 
is called the box dimension of U. 

The upper and the lower box dimension satisfy the following properties that can 
be looked up in e.g. [18] or [3] . The upper (lower) box dimension of the Cartesian 
product is the sum of the box dimensions of factors. The upper box dimension 
satisfies the finite stability property: the box dimension of the finite union is equal 
to the biggest of dimensions. Furthermore, they both satisfy monotonicity property: 
the dimension of a subset is smaller than or equal to the dimension of the whole 
set. The upper (lower) box dimension of an image of a set under a Lipschitz map 
is smaller than or equal to the box dimension of the original. In particular, box 
dimension is invariant to bilipschitz transformations. 

2. Main results 

In Theorem[2]we state the main result, the box dimension of a hyperbolic saddle 
loop. Theorenr[i] and Lemmas [l] and [5] are used in the proof, but they are also 
independent results. All proofs are in Section [3] 

Recall the asymptotics of the Poincare map on a transversal to the hyperbolic 
saddle loop not passing through the saddle from m , see 0 in Section |TJ The 
asymptotics by definition determines the box dimension of the orbit. Compare 
it to Theorem [lj In EL it was shown that the box dimension of an orbit of the 
Poincare map on a transversal away from the vertex does not reveal the codimension 
of the loop and the cyclicity uniquely. Here we compute the asymptotics on a 
transversal through the saddle. We see that, although the scale is slightly different, 
the information carried in the asymptotics is the same. 
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Theorem 1 (Poincare map on a transversal through the saddle). Let t be a 
transversal through the saddle of a hyperbolic saddle loop of codimension k £ N. 
Then the Poincare map s ha P(s) defined on r has the following asymptotics, as 
s —» 0: 


i s 7 ’, r > 1 the ratio of hyperbolicity, if k = 1, 
s k ~ 1 , if k even, 

s fc (—logs), if k odd, k > 1. 


Corollary 1. The box dimensions of orbits of the Poincare maps on a transversal 
away from the saddle vertex and through the vertex are equal to, respectively, 


dims 


1 

L k 1 


1 - 


2 

k+l ’ 


k even, 
k odd. 


dims = 



k even, 
k odd. 


Along the loop, we have two distinct directions with different rates of growth of 
Poincare maps: transversals through and not through the saddle. The asymptotics 
is slower through the saddle. This is new, compared to the focus and limit cycle case 
treated in m, and similar to degenerate and weak foci cases [3]. We check here if 
the saddle influences the box dimension of the trajectory, that is, if the accumulation 
of density of the trajectory lies at the saddle or away from the saddle. 

Let Xo be an initial point lying close to the loop and let S(x 0 ) denote the spiral 
trajectory with the initial point Xq, accumulating at the loop. 


Theorem 2 (Box dimension of the spiral trajectory around a saddle loop). Let 
k > 1 be the codimension of the saddle loop. Then 


dim B (S(x o)) = 


2 -|, 


2 - 


2 

fc +1 > 


k even, 
k odd. 


Note that the dimension of a trajectory from Theorem[2]is the product dimension 
around any transversal not passing through the saddle (by Flow-box theorem, see 
Section [l]) . Given the finite stability property of box dimension, the accumulation 
at the saddle is obviously not dense enough to carry the box dimension. 


Sketch of the proof. We compute the box dimension dividing the trajectory in 
two parts, S(xq) = Si(xq) U S%( xq), and using the finite stability property of box 
dimension, Sj^xo) is the non-regular part of the spiral trajectory near the saddle, 
between the transversals {x = 1} and {y = 1} (a family of hyperbolas). <S2(a;o) is 
the remaining regular part of the trajectory. We first compute the box dimension 
of Si(xo). The main tools are Lemma [I] and Lemma [2] stated below. First, in 
Lemma [l}_we show that Si(xo) car b by a bilipschitz mapping, be transformed to 
a simpleij_j family of hyperbolas, intersecting the transversals at the points with 
the same asymptotics as the original family of hyperbolas. Then, in Lemma [2j we 
compute the box dimension of a family of hyperbolas parametrized by a discrete set 
with the known box dimension on a transversal. To compute the box dimension of 
the regular part S^xo), we can directly apply the Flow-box theorem from Section[lj 
In detail, the proof is given in Section [3| 

The statement of the following Lemma[l]is the analogon of the flow-box theorem, 
but for a flow near a singular point whose portrait is a family of hyperbolas. It 
concerns straightening the original family to a simpler family of hyperbolas. 


'Lee Definition |Tj 
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Lemma 1 (Orbital linearization). Let X be an analytic vector field with a saddle 
at the origin: 


(3) 


x = x + P(x,y), 
y = —r • y + Q{x,y ), r > 0. 


Then there exists a neighborhood of the origin and a local diffeomorphism acting 
quadrant-wise, that transforms the phase portrait of ([3]) to the phase portrait of its 
linear part. 


Recall the C r -normal forms for hyperbolic saddles stated in [151 Theorem 13]. 
The field ^ above is not necessarily CMinearizable in the sense of C 1 -conjugacy 
of fields, but its phase portrait is, by Lemma[l] (^-equivalent to the phase portrait 
of its linear part. It is a weaker statement than C 1 -conjugacy of fields, see e.g. [T5] 
or m for precise definitions. 


Definition 1 (Family Tl r ,s )• Let r > 0. By 

LL r ,s = {(*, y ) G (0,1] x (0,1] | x r y = c, c G S}, ScM, 
we denote a family of hyperbolas^ parametrized by S Cl. 

Note: The family TL r ,s i n Definition [I] consists of integral curves of a linear field 
in the first quadrant: 

{ x = x, 
y = —r ■ y. 

Note that the set S and S l ^ r = {c 1//r | c G 5} respectively consist of points on the 
transversal {x = 1} and {y = 1} respectively, where the family of hyperbolas TL r ,s 
intersects those transversals. 


Lemma 2 (Box dimension of TL r ,s)‘ Letts' C (0,1] be a sequence of points ac¬ 
cumulating at the origin, with [^] distances between the points eventually decreasing. 
The box dimension of the family of hyperbolas 

LLr.s = {(ar, y) G (0,1] x (0,1] | x r y = c, c G S} 


is equal to 


dim B (LL r ,s) 


1 + dim b S, r < 1, 

1 + dims(S' 1//r ), r > 1. 


Note that, by Lemma[2j the box dimension of the family H r $ is in fact the bigger 
of the two box dimensions of parallel flows at transversals {x = 1} and {y = 1}. 
The accumulation at the saddle point is not what prevails in box dimension. The 
accumulation of density is not at the vertex. 

We illustrate the product statement of Lemma[2]on Figure[l]below. In the figure, 
r > 1. The family of hyperbolas {x r y = c, c G S}, intersects the transversals 
{x = 1} and {y = 1} in one-dimensional sequences S and S 1 / r respectively. Box 


hi set of level curves of function F(x,y) = x r y 

®Every mention of the interval (0,1] can be replaced by interval (0, <5], S > 0. The analysis is 
local, in a neighborhood of the singular point (the origin). Here, (0,1] X (0,1] is taken just for 
simplicity. 

^except at most finitely many first distances 
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dimension of the sequence S 1 / r is bigger than the box dimension of the sequence S 
(, x r for r > 1 is a Lipschitz map). The accumulation of density is therefore around 
the horizontal transversal {y = 1}, and the set of hyperbolas takes the product box 
dimension around {y = 1}, that is, 

dim s ('H r ,s) = dini B (S' 1/V ) + 1. 



FIGURE 1. Family T-i r ,s of hyperbolas in Lemma[2] r > 1. 


In Example [l] we give examples of sequences S generated by Poincare maps on 
transversals of hyperbolic saddle loops. 


Example 1 . Let g : (0,5] —>• (0,5], 5 > 0, and let f = id — g. For 0 < so < 5, we 
consider the orbit of g with initial point Sq : 


S = {s„ = </ on (s 0 )| n e N 0 } C (0,5]. 


In cases (1) and (2) below, the fixed point s = 0 is parabolic, and in case (3) it is 
weakly or strongly hyperbolic. A saddle loop with the ratio of hyperbolicity of the 
saddle r > 1 has Poincare maps of strongly hyperbolic type. If the ratio is r = 1, the 
Poincare maps are either of parabolic or of weakly hyperbolic type, see H5I Chapter 
5], 

(1) Let 

f(s) ~ s a , s —> 0; a > 1. 

By m Theorem 1], 


_ 1 

s n ~ n a ~ 1 , n —y oo, 


dim# 5=1-. 

a 


( 2 ) 


Let 

f(s) 

By P3J Theorem 2], 


s Q (—logs), s —> 0; a > 1. 

dims S = 1-. 

a 
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(3) Let 

g(s) = ks + o(s), 0 < k < 1 , or 

g{s)=Cs p + o(s p ) ) 0>1, C>0, s-> 0 . 

The orbit S accumulates at zero exponentially fast. There exist 7 £ (0, 1) 
and C > 0 such that 

0 < * n < cy\ 

By® Lemma 1, Theorem 5], 

dims S = 0. 


3. Proofs of main results 


Proof of Theorem [7} We mimic the proof of Poincare asymptotics on a transversal 
not passing through the saddle from (T5J Section 5], adapted for transversals through 
the saddle. Since the asymptotics is invariant to conjugacy, it suffices to work with 
finite C fc -normal forms of the vector field from [Tol Theorem 13], k £ N. 

Let r denote a transversal through the saddle. For simplicity, we take the diago¬ 
nal t = {y = x }. By Di and D 2 , we denote the transition maps from the horizontal 
transversal {y = 1} to r and from r to the vertical transversal {x = 1}. Then, 
D = D 2 o Di is the saddle transition map. Let R denote the regular loop transition. 
The Poincare map on r is given by the following composition: 


(4) 


P = Di o R o D 2 . 

Let r = 1. By m , after change of variables to (x,u), u = xy , 

\x = x, 


E oo 


From Q it can be computed that D(s) = u(— Ins, s), see [To] Section 5.1] for 
details. Similarly, we compute: 


(5) 


D 2 (s) = u{— In -^=, y), 

sV 2 


By m Theorem 14], 

u(u, s) = s + as k LO + h.o.t., a £ K, k £ N, k>2. 


Here, w is a logarithmic monomial and the higher order terms are meant with 
respect to a lexicographic order imposed on monomials s m Lo n , m, n £ N. Hence, 

D(s) = s + as k (— logs) + h.o.t. 

Suppose that R(s) = ais + a 2 s l + o(s l ) is of order l, l > 2. The Poincare map 
Pi = Ro D on the transversal {y = 1} obviously has the asymptotics: 

( 6 ) Pi(s) = ais -+[ 2 ]min{a 2 S i , aais fc (— Ins)} + h.o.t. 


7 


Here, min is meant in the sense: of the lower order with respect to lexicographic order. 
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Furthermore, by 0, we get: 

D 2 (s) = y + y-s 2fc (— In s) + h.o.t., 

Di{s) = V2S 1 * 2 H- < ^-=s k ~ 1 ^ 2 (— Ins) + h.o.t. 

2 y 2 

Componing P = Di o R o P> 2 , we get that 

(7) P(s) = ais + min{a 2 ai2~* • s 21 ^ 1 , aa\2~ k • s 2fc_1 (— Ins)} + h.o.t. 

Compare § and 0 , and the statement of the lemma for r = 1 follows. 

Now let r > 1. If the saddle is nonresonant (r ^ Q), it is by [T5I Theorem 13] 
locally C^-linearizable, k £ N. The linear system 

{ x = x, 
y = ~ry 

has an explicit solution x(t) = Ce 4 , y(t) = De~ rt . Using explicit formulas, we 
easily get that 

D(s) ~ s r , Di (s) ~ s^T 1 , D 2 (s) ~ s r+1 , s —> 0. 

Therefore, Pi(s) ~ s r , P(s) ~ s r , s —» 0. 

On the other hand, let the saddle be resonant (r > 1, r £ Q). We put r = p/q. 
Similarly as in the first case, after the change of coordinates to (x, it), it = x p y q , we 
compute 

u(cu, s) = s + as k co + h.o.t., k € N, fc > 2 . 

Since D(s) q = u(— lns,s p ) = s p + o(s p ), it follows that D(s) ~ s r , s —► 0. We 
compute 

£) 2 (s) = u( - In j ) /? ~ s r+1 , s -A 0 . 

From D = D 2 oDj, we have that .Di(s) ~ s?+i, s —>• 0. Thus, 

Pi(s)~s r , P(s)~s r , s->0. 

□ 

Proof of Corollary [Tj It follows directly from Theorem |T] and Example [l] □ 

Proof of Lemma[ 7] Let X have a nonresonant saddle at the origin, r Q. By [T51 
Theorem 13], such a vector field is linearizable in a neighborhood of the origin: it 
is even diffeomorphically conjugated to its linear part. This proves the statement 
in the nonresonant case. 

We prove here the more complicated, resonant case. Although not necessarily 
linearizable, we show that the phase portrait is quadrant-wise diffeomorphically 
equivalent to the phase portrait of the linear part. Let X have a resonant saddle 
at the origin, with r € Q+. Let r = | with p, q £ N, (p,q) = 1. By Theorem 
13 in T5], there exists an integer N G N, the coefficients a 2 , ■ ■ ■ ,on+i € R, and 
a neighborhood of the origin, such that X is diffeomorphically conjugated to the 
finite polynomial vector field 

x = x, 

= -r-y + l Eh «i+i • {x p y q Y ■ y- 


( 8 ) 


y 
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We construct a diffeomorphism F(x,y), acting quadrant-wise in a neighborhood 
of the origin, which sends phase curves of field ([8]) to phase curves of its linear part. 
We show here the construction of F 1 (x, y) in the first quadrant. Afterwards we glue 
functions F I,II,III,IV constructed in each quadrant to a global diffeomorphism F 
at the origin. 

We proceed as in [15, 5.1.2] (a similar technique was used there for obtaining 
the Dulac map at the resonant saddle). We solve the system ([8]) by substitution 
u = x p y q , whereas we get the system 


(9) 


x = x, 

u = Li =2 ab ¬ 


solving by expanding u(t,Uo) in series with respect to the initial condition 
uq, we get that 


( 10 ) 


N 

u(t, u 0 ) = «o + E 9i(t) u b- 

i—2 


The form of g?;(f), i > 2, is described in Proposition 10 in [l~5’l : gi{t) are polynomials 
in t, of degree at most i — 1. Therefore, we easily obtain the bounds: 

(11) \gi(t)\ < Qtf- 1 , |ff'(t)| < Af“ 2 ; Q, Di > 0, i = 2,..., N, 


for t sufficiently big. 

Let T\ = {y = 1} be a horizontal transversal to the saddle. We should in fact 
consider the transversal at some small height <5 > 0 instead at height 1, but the 
computations are the same. Using (10) and u = x p y q , we can now derive the 
formula for the phase curve of the cut-off field ( |8|), passing through the initial point 
(ar(0), 2/(0)) = (s, 1) £ t\ . We put uq = s p in (|10|) and solve x = x. We get that 
t = log -, and then, for the phase curve through (s, 1), we have the formula: 


( 12 ) 


S / 

V = w( 1 + 


N 

E^g-)* 

z — / s 

2=2 


.p(i-l)\ 1 /9 


S < X < 1 . 


5 _ 


The phase curve of the linear part passing through (s, 1) is, on the other hand, 
given by 


(13) 


y=—, s<x< l. 
x r 


We now define the mapping F 1 of the first quadrant (0,1] x (0,1] to itself, sending 
phase curves of ([8]) to phase curves of the linear part in the following manner. For 
any point (x, y) £ (0,1] x (0,1] close to the saddle, there exists a unique phase 
curve of the linear field passing through it, and it is determined by the point (s, 1) 
on T\. We consider the phase curve of ^ passing through the same point (s, 1). 
F I {x,y) is then defined as the orthogonal projection of (x,y) on this phase curve. 
More precisely, by (12) and we get the formula for F I (x,y): 


N 

(14) F I (x,y) = l X , y (i + E^( 1 °g(i/“ 1/r ))-^ (i “ 1) -y 9(i " 1) ) 


1 / 9 ' 


i=2 


(x,y) £ (0,1] x (0,1]. 


Function F 1 is obviously well-defined, continuous and differentiable in (0,1] x (0,1]. 
We can obtain similar formulas for F II,III,IV in other quadrants. In the first and 
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in the second quadrant, in (13), we have y = s r /x r . In the third and the fourth 
quadrant, we have y = — s r /x r . By glueing the quadrants together, we get function 
F, defined on the unit square without the coordinate axes, by the formula 


F(x,y) = 


(x, y(l + 


N 

E 

i=2 


gi{\og{\y\- x / r ))x^-^y^-^ 


1/9 


(15) (x, y) G ([-1,1] X [-1,1]) \ ({x = 0} U {y = 0}). 

It can be checked that the functions can be extended to the coordinate axes to a 
continuously differentiable function F(x,y) defined on the square [—1,1] x [—1,1], 
Then we apply the inverse function theorem at the origin. We conclude that F is 
a local diffeomorphism. Moreover, by construction, it leaves the axes invariant and 
maps each quadrant to itself. 


Let us extend F from (15) to the coordinate axes, and check that the obtained 


function is continuously differentiable in [—1,1] x [—1,1], We define 

F{x, 0) := (x, 0), x G [-1,1], F( 0, y) := (0, y), y G [-1,1]. 


Note that F cannot be extended continuously simply by formula (14) to y = 0 


due to the logarithmic term. However, using (11), lim^o F(x, y) = 0, x G R. 
Therefore, extended as above, F is continuous on [—1,1] x [—1,1], 

Furthermore, F given by (15) is differentiable on ([—1,1] x [—1,1]) \ ({& = 
0} U {y = 0}). We can show, by direct computation of the differential and using 
bounds (fTTj), that: 


(16) 


DF(x, y) = 


1 0 

d x F 2 {x,y) d y F 2 (x,y) 


, x^0, y ± 0, 


Here, 


lim d y F 2 (x,y ) = 1, lim d y F 2 (x,y) = 1, 

y —¥ 0 x —^0 

lim d x F 2 (x , y) = 0, lim d x F 2 (x, y) = G(y). 

■ y —> 0 x —>-0 


G (y) = («' y9+1 ' 92 ( l ° g \ y \ P=l ' 

V 10, p> 1. 


Obviously, G(y) —> 0 as y —Y 0. We can check directly by definition of differen¬ 
tiability at (x,0) and (0 ,y), using bounds (11) and (16), that F extended to the 


axes in the above manner is continuously differentiable at the axes and that the 
differential is given by 


(16), 


DF(x,y) = < 


1 0 
0 1 


(x,y) G ([-1,1] x [-1,1]) 

\({x = 0} U {y = 0}), 

x G [-1,1], y = 0, 


x = 0, y G [-1,1]. 


1 0 
G(y) 1 

Note in the course of the proof that C 1 is the best class that we can obtain applying 


bounds (11). 
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Proof of Lemma\^ The idea is similar to the proof of the box dimension of a Carte¬ 
sian product in the plane. There, the area of the e-neighborhood of a 2-dimensional 
product U x [0,1] can be expressed using the length of the e-neighborhood of the 
1-dimensional set U. 

Let r < 1. The family T~L r ,s is a family of phase curves (in the first quadrant) of 
the linear saddle with ratio of hyperbolicity r. By Flow-box theorem from Section [l] 
we get that the box dimension of P r ,s on a small rectangle around {x = 1 } is equal 
to 1 + dim B S. Around {y = 1 }, the box dimension is equal to 1 + dim s (S' 1//T '). 
Since y = x 1 ^, r < 1, is a Lipschitz map, we have that dirng(S ,1 / r ) < dim^ S. By 
the monotonicity property of box dimension, we get 

(17) dim b (Pt,s), dim B (P r ,s) > 1 + dim s S. 

If r = 1, the hyperbolas are symmetric with respect to the diagonal y = x. 
If r < 1, we first symmetrize the family by the change u = x r , v = y, which 
is inverse-Lipschitz. In the new coordinate system ( u,v ), we get the symmetric 
family of hyperbolas Pi t s = \ uv = c| c £ S}. We compute the box dimension of 
the symmetrized family Pi t s, and conclude that 

dim B P rt s < dim B P\,s- 

The box dimension of His is, by symmetry and by finite stability property of 
box dimension, equal to the box dimension of its subset between the diagonal y = x 
and the transversal {x = 1 }. 

We now estimate the area of the e-neighborhood of Pi t s in this area, estab¬ 
lishing the almost product relation with the length of the e-neighborhood of one¬ 
dimensional set S on the transversal {x = 1}. The intersection points of Pis and 
the diagonal {y = x} form the set y/2S 1 ^ 2 on the diagonal. The distances of the 
intersecting points on the transversal {cc = 1 } are smaller than the distances of 
corresponding points on any other transversal in the area, including the diagonal 
(at least close to 0 , considered on ( 0 ,5] x ( 0 ,5], S > 0 ). 

We compute the area dividing the e-neighborhood into tail and nucleus. By tail 
of the e-neighborhood, T 2 , we mean the disjoint neighborhoods of the first finitely 
many hyperbolas. The remainder of the e-neighborhood, where the neighborhoods 
of hyperbolas start overlapping, is called the nucleus , and denoted iV 2 . For the 
idea of division in tail and nucleus, see [18]. Since the distances are the smallest 
on transversal {a; = 1 }, the critical index separating the tail and the nucleus is the 
same as for one-dimensional set S on {x = 1}. Let Tf, N E denote the tail and the 
nucleus (in dimension one) of its e-neighborhood. We now bring Tf and T 2 and 
Nf and TV 2 into direct relation, and express the box dimension of P\,s by the box 
dimension of S. 

Since the lengths of hyperbolas are bounded from above and below in (0, <5] x (0, d], 
we get that|^] 

(18) A(T 2 ) ^\Tl\+e- \Tl\ ~ IT, 1 ), e -a 0. 


®The notation ~ stands for: / ~ g, as x —> 0, if there exist A, B > 0 and 5 > 0, such that 

f (x') 

A < < B, 0 <x <6. 

9W 
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By definition of box dimension, it holds that 

ra = 0 , lim I"* 1 ' 


(19) 


g-l — dim B S-S 


By ([18]) and ( |19[ ), we get 
A(Tf) 


lim 


= lim 


e-j-O £-1—dims S—S 

A(Tf) 


= 0 , for all 5 > 0 . 


= 0 , for all S > 0 . 


(^0) £ l-dim B S-S “4J £ 2-(dim fl S+lH 

For the nucleus, we give an upper bound on the area of the e-neighborhood. The 
hyperbola separating the tail and the nucleus has the equation uv = | TVj |. The 
area A(Nf) is smaller than or equal to the area between the tt-axis, the hyperbola 
H ne , the diagonal y = x and the transversal {x = 1}: 

A(Nf) < C\N}\ + De + f ® dy, e < e 0 . 

d l Mi i 1 / 2 y 

Take any fixed Jo > 0. There exists a small v > 0, such that <5 0 — v > 0. Integrating, 
it holds that there exist C > 0 and £o, such that 

A(N e 2 ) ^ C\Nl\(-loge) ^ C\N^ 


■e v (- loge), £<e 0 . 


£-1—dirn^ S—8o — £-1—dime S—8q — £-1—dim# S— (<5o —^ 

Passing to limit as e —> 0 in the above inequality and using (19), we get that 

A( N?) 


( 21 ) 


Am 

e— f 0 £l-dim B S-S E^0 g2—(s+1)—5 


lim 


= lim 


= 0 , for all S > 0 . 


By (20) and (21), it follows that dim B ('H lj 5 ), dim^^ys) < 1 + s for the sym¬ 
metrized family. 

It follows that dims('H r ,s) < dims('H 1 ,s) < 1 + s. The same for the upper box 
dimension. Using lfl7|)7 we finally get that, for r > 1, 

dim B (Hr, s) = 1 + s. 

Now let r > 1. Then, 1/r < 1. We have 

?U,s ={(x, y) £ (0,1] x (0,1] | x r y = c, c £ S} = 

={{x, y) £ (0,1] x (0,1] | y 1/r x = c, c £ S 1/r j. 

Changing the roles of x and y and substituting S for S 1 ^ , by the first part of the 
proof we get that 

dim B (’H ri 5) = 1 + dim B (5 ,1 / r ). 

□ 


Proof of Theorem Q Let the saddle loop be of codimension k, k £ N. Without 
loss of generality, suppose that the ratio of hyperbolicity is r > 1. The repelling 
case r < 1 can be considered as attracting by rescaling the time variable t —rt , 
whereas the phase portrait (and so the box dimension of a trajectory) remains the 
same. Take a spiral trajectory S(x o) accumulating at the loop, with xq close to 
the loop. We divide the trajectory in two parts: the part S) (x a ) passing the saddle 
vertex, lying in the square ( 0 , 1 ] x ( 0 , 1 ], and the remaining regular part S^x o) 
around the saddle connection. 

We first compute dims(S’i(a:o)). The asymptotics of the Poincare map s i —> 
P{s) on the transversal {x = 1} is given by ([!]). Let S C (0,1) denote the orbit 
of P with initial point So £ (0,1] - one of the points on {x = 1} where S(x o) 
intersects it. By Lemma [l] by a local diffeomorphism, the arcs of S±(x o) can be 
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transformed to a family = {(a;, y) £ (0,1] x (0,11 : x r y = c, c £ S'} of 
countably many hyperbolas. They intersect {x = 1} at FIS'. The box dimension of 
T~Lr,s is given in Lemma [2j using Corollary [Tj for the box dimension of S. Since a 
local diffeomorphism is a bilipschitz map, we have: 

dim B (Si(®o)) = dim?^ = 

It is left to compute the dimension of the remaining, regular part S^aio) of the 
trajectory. In this area, there are no singularities of the vector field. Therefore 
we directly apply the Flow-box theorem [ 8 ], see Section [l] and Corollary [I] The 
box dimension of S^^o) is computed as the box dimension of Cartesian product 
S x [0,4 l > 0: 


2 — k even, 

2 — ttzt, k odd. 


dim B (S 2 (so)) = 


i + (i-i) = 2- 


2 

k ’ 


2 r) = 2- 2 


k even, 
k odd. 


k+\> “ k+ 1’ 

Finally, by finite stability property of the box dimension, the result follows. 


4. Applications 

4.1. The cyclicity of a hyperbolic saddle loop. In recognizing the cyclicity 
of monodromic limit periodic sets for planar systems, one can use fractal analysis 
of its trajectories. By fractal analysis , we assume analysing the lengths of the e- 
neighborhoods of orbits of Poincare maps, as functions of e > 0. Also, analysing 
the areas of the e-neighborhoods of trajectories accumulating to the set. The box 
dimension of a set is, by its definition in Section [l] related to the asymptotic be¬ 
havior of the Lebesgue measure of the e-neighborhood of the set, as e —»• 0. The 
relation between the box dimension and cyclicity was first investigated in [19] , [20] 
and later in m- In weak focus and limit cycle cases (Poincare maps differentiable 
at 0 ) there exists a bijective correspondence between the box dimension of any 
orbit of the Poincare map and the cyclicity of a set. Further, the bijective corre¬ 
spondence was established between the cyclicity and the box dimension of a spiral 
trajectory around a focus or a limit cycle, using Flow-box theorem or its appropriate 
adaptation, see [ 201 . 

The paper m was motivated by the attempt to do the same for hyperbolic saddle 
polycycles, in particular for the simplest saddle loop. Here, the Poincare maps 
contain logarithmic monomials and are nondifferentiable at 0. Thus, comparison 
of the asymptotic behavior of the lengths of the e-neighborhoods to powers of e 
in the very definition of the box dimension is imprecise. See [13 p. 2502] or 
Example[l](l), (2) and note that the box dimension reveals the cyclicity ambigously. 

The asymptotic expansion of Poincare map on a transversal through the saddle 
vertex in Theorem [l] shows that, although slightly different, it carries exactly the 
same information as Poincare maps on any other transversal. Further, the box 
dimension of a two-dimensional spiral trajectory accumulating on the loop exhibits 
the same deficiency in revealing cyclicity, precisely formulated in Proposition [l] 
below. This is not surprising, due to the fact that Lemmas [T] and [2] develop a version 
of Flow-box theorem adapted for hyperbola-type flow, thus implicitely showing that 
all dimension information is already contained on a transversal. 


°By definition, only the density of accumulation of a trajectory on the loop is important for 
the box dimension of S(x o); the first finitely many windings can therefore be neglected. 



BOX DIMENSION, SADDLE LOOP 


15 


Proposition 1 (Cyclicity of a saddle loop and the box dimension of its spiral tra¬ 
jectory). Let (r, X\) be a genericp°| analytic unfolding of a (monodromic) hyperbolic 
saddle loop T. Let S(x o) be a spiral trajectory with initial point Xq accumulating at 
r, and 

dim s (5(x 0 )) = d£ [1,2). 

By the box dimension, the cyclicity of the loop is not uniquely determined. More 
precisely, 

CycHT.X,)e{^- d - 1 ,^}. 

Proof. By Theorem [5J a loop T with dimB(5'(a:o)) = d may be of codimension either 
2 PL 1 ~ 1 or Under genericity assumption on the unfolding, the codimension is 
equal to the cyclicity. □ 


4.2. Prospects: the box dimension of leaves of a foliation of complex 
resonant saddles. We consider germs of holomorphic vector fields in C 2 , with 
complex resonant saddles at the origin. That means, with the linear part of the 
form: 

jz = z, 

) w = —r ■ w, re Q+. 

Here, r e Q+ is the saddle hyperbolicity ratio. Put r = |, (jp,q) = 1. 

Such germs are either formally orbitally linearizable or formally orbitally equiv¬ 
alent to the germ 


( 22 ) 


z = z, 

w = w(—r + 


u k+1 \ 
1+Au fc >> 


u = z p w q , 


for some k e N, A e C. This result can be looked up in e.g. jl(n Sections 1,2,22], 
[2 Chapters 4,5], [T7]. 

Motivated by )13l , where the formal type of a germ of an analytic diffeomorphism 
is read from fractal properties of one orbit, our goal is to see if the box dimension 
of a leaf of a foliation can reveal the formal type of the saddle. 

Let L a denote a leaf of a foliation through a point a £ C 2 sufficiently close 
to the saddle. Let T\ = {w = rco} and r 2 = {z = Zq} denote a (two-dimensional) 
horizontal and a vertical cross-section. Let h w (z) denote the holonomy map induced 
by L a on t\ and h z (w) on r 2 . By im around each cross-section, a leaf of a foliation 
has a locally parallel structure (it is a bilipschitz image of a family of unit complex 
discs). 


genericity , we mean the following regularity condition on the unfolding from [6] or Qoj. 
Let Poincare maps for the unfolding (Xa) have the following asymptotic expansion, as s —> 0: 

P\(s) = ao(X)u 0 (s , A) 4-... + a k (X)u k (s, X) + h.o.t A G A C T. 

Suppose r is of codimension k < N, P(s, Ao) = ot k (Xo)u k (x, Ao) + h.o.t. That is, k conditions are 
imposed on the loop (c>!o(Ao) = . . . = a k _i(Xo) = 0). The regularity condition is that the matrix 

/ daj(\ p) \ 

V dXj ' i=Q,...,k— 1; j=l,...,k 

is of maximal rank k. That is, by the implicit function theorem, the coefficients ao(A),..., a^_i(A) 
can be freely chosen in X\. Under this assumption on the unfolding, the cyclicity of the loop is 
exactly equal to its codimension. Without this assumption, it could be smaller. 

■^Induced foliations are formally conjugated. 
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With respect to formal linearizability, there are two types of resonant saddles, 
see m Lemma 22.2]: 

(1) Formally orbitally linearizable saddles. It holds that 

h° q = id , h \= id. 


Orbits of holonomy maps on cross-sections consist of finitely many points. 
The box dimension of orbits is thus equal to 0. By product structure, we 
conclude that the box dimension of L a locally around each cross-section is 
equal to 2. 

(2) Formally orbitally nonlinearizable saddles. The iterates of the holonomy 
maps are complex parabolic germs, 

h° q (w ) = w + a\W kq + h.o.t., h° p (z) = z + b 1 z kp + h.o.t. 

Orbits of h z , h w , denoted S hz and S hw , consist of finitely many ( q , p) dis¬ 
joint orbits of h ° q , /i°f respectively. By m and by finite stability property 
of box dimension, 


dim B (S^ 2 ) = 1 — 


kq + 1 ’ 


dim s ( S hul ) = 1 — 


1 

kp + 1 


Here, k, p, q are as in (22). By product structure, we conclude that the 
box dimension of L a locally around cross-section n is equal to 3 — kq+i • 
Locally around 72 , it is equal to 3 — k f 1 . 


By monotonicity property of box dimension, under the assumption that box dimen¬ 
sion exists, we get: 


(23) 


dim s (L a ) > 


2 , 

max { 3 ^ T^+I’ 3 ^ ’ 


for case (1), 
for case (2). 


To verify the other side of the inequality (23), we need to compute the box di¬ 
mension of a leaf in a small neighborhood of the origin, where the product structure 
is lost. This is left for future research. Motivated by results of Lemma [2] for the 
planar saddle, we conjecture that, also in the complex case, we have equality in 

(231. 
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